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Abstract 

In the framework of the rate-independent large-strain Cosserat the- 
ory of plasticity we calculate analytically explicit solutions of a two- 
dimensional shear problem. We discuss two cases where the micro- 
rotations are stationary solutions of an Allen-Cahn equation. Thus, 
for a certain parameter range, patterning arises and the domain is 
partitioned into subsets with approximately constant rotations. This 
describes a possible mechanism for the formation of grains and sub- 
grains in deformed solids. 
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1 Introduction 



Experimental evidence suggests that the size of grains and subgrains in plas- 
tically deformed solids is not random, but rather relates inversely to the flow 

stress, 

a ~ zr 1 / 2 , 

see [3 E]- The foundations of this heuristic law are a long-standing open 
problem in materials science and solid mechanics. The qualitative and quan- 
titative prediction of D is particularly vexing in dynamic recrystallisation 
(DRX), where new grains essentially free of dislocations nucleate out of a 
highly fatigued sample. The experimental studies jH] on DRX show that the 
recrystallized grain sizes can be linked to the subgrain size although the de- 
pendency of dynamically recrystallized grain size and subgrain size on flow 
stress is different. 

Derby and Ashby [7], and Roberts et al., [18], addressed the grain size 
by making assumptions on the relevant criteria like back stresses at grain 
boundaries or imbalance of subgrain size across prior grain boundaries but 
did not touch the issue of mechanical behavior. A survey on the effect of 
diffusion on grain boundaries can be found in pp. 

Evidently, despite extensive research during the past decades, the physical 
mechanisms of grain formation are not yet understood. The subject is not 
only of academic interest but is also of immense importance for industrial 
processing, see e.g. [21 [101 [3J [2TJ, since treatments like hot rolling, cold 
rolling and annealing require reliable predictions of the microstructure and 
the texture evolution. 

In this article, without spurious assumptions or introducing an artificial 
mechanism, we show in a mathematically rigorous way why patterning in de- 
formed solids may occur. To develop our argument we compute analytically 
in special cases the state of a solid that is deformed along its boundary. We 
formulate the problem within the finite-strain Cosserat model of plasticity. 
The interesting feature of this model, unlike other established approaches in 
visco-plasticity, see the research articles and surveys [TjJJ [2UJ EH EE2] , is that 
it is a gradient model, i.e. a length scale is introduced. The results given here 
substantiate this property: it is demonstrated how for a certain parameter 
range the analytic solutions naturally lead to transition layers between areas 
of approximately constant micro-rotations R c , a parameter that specifies the 
local orientation of the material. 
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This paper is organised in the following way. In Section [5] we revise the 
rate-independent finite-strain Cosserat model as needed later and formulate 
an energy minimisation procedure that allows to obtain the time-discrete 
solution. We then restrict the plastic deformations to a priori given single- 
slip systems and adopt the formulation accordingly. In Section HI we apply 
the resulting model to a simple shear problem in 2D. The two-dimensional 
setting is chosen as this permits the representation of the rotations by one 
scalar parameter. In the following section [5] some exact solutions to this shear 
problem are computed with focus on the properties of the micro-rotations 
R c . As main result, two non-trivial examples are discussed where R c is 
connected with the solutions of a stationary Allen-Cahn equation. We end 
with a discussion of the results. 

2 The rate-independent finite-strain Cosserat 
model of visco-plasticity 

Let fl C M. d be a bounded domain with Lipschitz boundary that serves as a 
reference configuration. The deformation of the material is controlled by a 
mapping <p(-,t) that maps fl diffeomorphically to the deformed state fit at 
time t > 0. Because of <p(-, 0) = Id it holds det(D<p(t)) > for all t > 0. 

The gist of the Cosserat approach is to multiplicatively decompose the 
deformation tensor F := Dtp into a plastic part F p and an elastic part F e by 
virtue of 

F = F e F p , (1) 

and to split F e by 

F e = R e U e (2) 

into a stretching component U c G GL(R d ) and a rotation part R e G SO(d), 
where 

SO(d) := {R G GL(R d ) | det(R) = 1, R l R = Id} 

denotes the special orthogonal group. Equation (JTJ) states that the plastic 
deformation of the material precedes the elastic deformation. In general, 
U e need not be symmetric and positive definite, i.e. (|2J) is not the polar 
decomposition. 

Introducing the dislocation density k, frame indifference (which amounts 
here to the invariance of the energy under rigid body motions) implies that 
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the mechanical energy density is the sum of three functionals, 



W(F e , k) = W st (U c ) + W c (K e ) + V(k), (3) 

with the stretching part W st , the curvature part W c , and the energy due to 
immobilised dislocations V. With K c := R l c D x R c = (_R*9a;;i?e)i</<d we des- 
ignate the third-order (right) curvature tensor. Eqn. ([3]) tacitly assumes that 
the mechanical energy depends on the elastic energy F D only. For arbitrary 
materials, this need not be the case. 

Following the ideas in [TB], we now derive a time-discrete formulation that 
allows to compute the evolution of the material by minimising the mechanical 
energy. For known values (F po , k,q) of the previous time step and discrete time 
step h > 0, the values of (<p, R c , F p , k) at time t + h are calculated. 

With P := F^ 1 , Pq := F^ 1 we choose the discretisations 

Ahtj? \ ._ Id ~ P ~ lp o a h K-Kp 

d t {F v ).- , d^.-—^. 

By fcxt(t) we denote external volume forces applied to the crystal body, 
^ext(^) are external volume couples. Let T D C dQ be the Dirichlet bound- 
ary (for simplicity postulated invariant in time) which we assume to be a 
smooth curve with positive (d — l)-dimensional Hausdorff measure. On 
we prescribe the Dirichlet boundary conditions 

Wd = 9d, R c \r D = Rd- (4) 

Since we are only interested in the traction-free case, we henceforth set := 
dVt. The time-discrete mechanical energy functional for given (F Po , kq) then 
reads 

E(ip,R e: F p: K)(t) := J [w st (U c ) + W C (K C ) + V(k) - f cxt (t)- V - M cxt (t):R e 



hQ*(d h t (F p ) : dlK) dx. (5) 



The term /iQ*(d t h (F p ), d£n) in © represents the dissipated mechanical energy 
in a discrete time interval of length h, defined by the Legendre-Fenchel dual 

Q*(F p ,k):= sup{X:F p + 7iK-Q(X,n)} (6) 
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of the plastic potential 



0, for Y(X, tt) < 0, 



and y is the yield function 

Y(a, it) := ||devsym<7|| — ay — n < 0. (7) 

This is the von Mises approach, devu := a — 4tr(er)Id designates the de- 
viatoric part of a tensor a, symcx := 0.5(er + a*), and ay > is the yield 
stress. 

The solution of the minimisation problem 

E(tp, R e , F p , k) -> min (8) 

subject to the boundary conditions (jlj) yields (<p, R c , F p , k) at time t + h. 
The first variation of E w.r.t. the (implicit dual) variables (X, tt) leads to 
the time-discrete plastic flow rule (see [19] for details) 

(F po ^(F p ),^)e9 sub g(X,vr). (9) 

The variation of E w.r.t. (tp, R e ) gives back the field equations. The variation 
w.r.t. (F p , k) returns the definition of the back stress X and the hardening 
modulus 7r, 

X = ^— , n=-V(K). 



3 Single-slip systems 

We specialise the minimisation problem (JSJ), (JSJ) to single-slip plasticity. In 
[1], it had been shown that the relaxed behavior of crystals restricted to 
single-slip under the infinite latent hardening constraint is identical to the 
dynamics of materials with general multislip plasticity. 

Let m a G M d be given slip vectors, n a G M d be given slip normals with 
Irrio, | = \n a \ = 1, m a -n a = for 1 < a < I, I the number of slip systems 
present. Let 

Y(a, tt) := max < \m a -an a \ > — ay — ix 

l<a<I L J 
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be the corresponding yield function. For single-slip systems, 



F p = F p (j) := Id + ^ lam a ®n a , 



0=1 



with the parametrisation 7 = (7 a )i< a <r £ ^ ■ 

As can be checked, the dissipated energy satisfies the identity 



Q*(A,k) 



°V ELl if ^4 = I]a=l7a m a® ra a and So=l lTa| + & < 0, 



oo, 

Therefore (18]) becomes 



else. 



£(</?, i? c , 7, k) 



W Bt (RlDlfFp (7) _1 ) + W C (K C ) + - /exf¥? 



-M cxt : i? c + a Y 1 7a - 7°l dx min ( 10 ) 



a=l 



subject to /i(7, «) := Y^l=i 1 7a - 7° I + K _ K o < and 
From now on, we set M ext = 0, f ext = 0. 



4 The shear problem in 2D 

In general, even in two space dimensions the minimisation problem ( TTUj) 
is too complex to be solved analytically and numerical computations are 
required. In order to proceed, we make two simplifying assumptions. The 
first assumption is that there is only one active slip system, 1 = 1, such that 

F p (7) = Id + 7 m®n, P(l) — Id — 7 m®n. (11) 

The analysis can be carried out without this assumption, but the formulas 
become very lengthy. Secondly, we assume that the total deformation is a 
shear of the material along this slip system, 

Dtp(t) = Id + /3(t)m®n in H (12) 

for a given scalar function 0(t). On dQ, this identity must be imposed by 
choosing in (J3J) the appropriate boundary conditions on if. Eqn. (TT2"j) states 
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that the deformation at every point in f2 follows this prescribed deformation 
at dQ, i.e. postulates the validity of the Cauchy-Born rule. 

Assuming (|12p . the mapping 7 1— > W s t(R l c D(pP (7)) is convex (see Lemma[T] 
below). This implies that the constraint /i(7, k) < is satisfied with equality, 

|7 - 7o| + « - = 0. 



Therefore we can resolve this condition and (110]) can be written as the un- 
constrained minimisation problem 

^(i? c , 7 ) = [ rty st ( J R*(Id+(/3- 7 )m®n)) + W c (ii- c ) + y(«: -|7-7 |) 



dx — > min, (13) 



n 

+a Y \l - 7° I 

Re\an = Rd- 



To explicitly compute the solutions to (fT3|) . we pick 

V{k) := qk 2 . (14) 
For W st , W c we make the general ansatz, cf. [T5j . 



W st (U e ) := /i||symt/ c -Id|| 2 + /i c ||skwt/ c || 2 + ^|tr(f/ c -Id)| 2 , (15) 

W C (K C ) := ^^(l + a^HA'ef) 

x(a 5 ||symA c || 2 + a 6 ||skwA c || 2 + « 7 |tr(i^ c ) | 2 ) ^ , (16) 

for positive material parameters g, /1, fi2, A, /i c , constants L c > 0, «4 > 0, 
«5 > 0, «6, OL7 > 0, p > 0, q > 0, and where skwA := 0.5(A — A f ), 

\\A\\ := y/tr(A t A) 

is the Frobenius matrix norm, and tr(A) := J7 An the trace operator. Set- 
ting L c := a/2, a 5 = a 6 := 1, a 4 = a 7 := 0, p := 1, the definition ffT6]) 
simplifies to 

W C (A C ) =/i 2 ||A e || 2 . 
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Since 



\K e \\ 2 = ^W^IUW 2 = ^2tT{d Xl H^R e Hid Xl R e 
1=1 1=1 
d d 

= 5>(^j£0 BI jy = x;ilS.,fl- 



|2 



Z=l 



the definition of W c simplifies further to 

W C (R C ):= fi 2 \\VR c \\ 2 . (17) 
In two space dimensions, we can benefit from the explicit representation 



R e = R e [a) = 
Finally, by direct computations 



cos a — sin a 
sin a cos a 



a E [0, 2vr) 



;is) 



ll^el 

||sym{7 c -Id| 



sm a — cos a 



d Xl a 



cos a — sm a 
||sym[ J R c (-a)(Id + (/3-7)m<gm)] - Id|| 2 
2(1 — cos a) 2 + 2(/3— 7)(1 — cos a) sin a 

+ (/3— 7) 2 ^ sin 2 a + - cos 2 , 



2\d Xl a\ 2 , Z = 1,2, 



|skwf/ c || 2 



2 sin 2 a + 2(/3—y) sin a cos a + (/3—y) 



2 cos a 



|tr(j7 c -Id)| 2 = 4(l-cosa) 2 + 4(/3-7)(l-cosa) sina + (/3~7) 2 sin 2 a. 

In summary, using the definitions ( fT4l) , (fl~5l) , ( JT7l) and letting Ko = 7o = for 
simplicity, ( 113]) becomes (with /^(od) = i?£>) 

Ep(a,y) = J 2/i 2 |Va| 2 + f?7 2 + <Ty|7| + 2(A+/i)(l — cosa) 2 + 2/i c sin 2 a 
n 

— 2(7 — / 9)((/i c — A — /i) cosa + A+yu) sina 



+ — 7 (// + (A + /i) sin 2 a + [i c cos 2 a) dx — > min, (19) 



ot\m = ot D . 



(20) 
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Lemma 1 The mappings 7 1— > Ep(a,~f) and 7 1— > W s t(U e (a, 7)) are convex. 
Proof From (1191) we have 
d 2 W st (U e (a^)) 



«9 7 2 



/i + (A + /i) sin a + /i c cos a > 0. 



Since V(7) and 7 (-)■ cry|7| are convex, this implies the convexity of Ep(a, ■). 

□ 

We study (HSJ, for a e ff 1 ' 2 ^; [0,2tt)) and /3(t), 7 e L 2 (fi; R), 
where f3(t) is given and controls the shear of Q. From the direct method in the 
calculus of variations, since E is coercive and weakly lower semicontinuous, 
the existence of solutions (a, 7) follows. 

5 Analytic solutions of the 2D shear problem 

In this section we study properties of the solution (0,7) of ( fT9l) for several 
choices of the parameters and are especially interested in cases where tran- 
sition layers in R e occur. 

5.1 The limiting case fi2 — > 00 

This limiting case corresponds to rigidity in bending. After rescaling, only 
the curvature energy W c contributes and (fT9l is equivalent to 

E(a) = 2 / I Val 2 dx ->■ min 



n 



subject to ()20p . The minimiser a solves the classical harmonic problem 

—Aa = in Q, a\sn = an- 

5.2 The limiting case g, ay \ 

The limit g \ 0, oy \ corresponds to an ultra-soft material where the 
dislocations do not contribute to the mechanical energy. Here, Equation ( |T9|) 
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simplifies to 

Ep(a,y) = y 2/i 2 |Va| 2 + 2(A+/u)(l-cosa) 2 + 2;U c sin 2 a 



-2(7 — ( 5)((/x c — A— /x) cosa + A+/i) sina 
(7-/3) 2 



-(/x+(A+/i) sin a+/i c cos a)) dx — > min .(21) 



Clearly, due to convexity in 7, the optimal 7 follows the given deformation 

7o P t(0 = 

Computing the Euler-Lagrange equation of (|2Ti) . we need to find a bounded 
solution a opt G [0, 27r) of the semi-linear elliptic partial differential equation 

— fi2 Aa + [A + fi + (/i c — A — fi) cos a)] sin a = in Q (22) 

with boundary conditions fl20l) . This equation has been studied before in 
[T7] . Based on a separation ansatz, special solutions have been found in [T3] , 
First consider the unphysical case /i2 = 0. Then evidently 

a op t e {0,7r} in Q, a opt \ dn = a D . (23) 

Since there is no regularising term |Va| 2 in E, a opt G {0,7r} a.e. in Q and 
a G L 2 (Q; [0, 2tt)) has no additional regularity 

When fi2 > 0, the solutions to (1221) are continuous and stationary (in 
time) solutions of the Allen-Cahn equation 

d t a = /i 2 Aa - J' (a) (24) 

with the nonlinearity 

J(a) := — (A + ji) cos a + — — — sin 2 a. (25) 

J has a local minimum at 0. For 

/i c >2(A + /i), (26) 

there is exactly one other local minimum at a = n as illustrated in Fig. [TJ 

The properties of the solution to (l24"j) when fl26l) holds are well known: a 
tends pointwise to one of the minimisers of J (here and 7r), but transition 
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Figure 1: Plots of J for A = pt, — 1. Left: /i c = 12. J has the two local 
minima and it. Right: fi c = 1. J has only one minimum at 0. 



layers of width yfpv and slope C/^/JI^ occur. For small /i2, the interface 
moves by its negative mean curvature. For t — > oo, a becomes stationary 
and when periodic boundary conditions are applied, a = or a = it in Q. 
If Dirichlet boundary conditions (|2"U1) are imposed, a sufficient condition for 
the presence of transition layers in the stationary limit t — > oo is 

/ 7T\ /3lT \ 7T 3lT , . 

[0<a D <-j or ^Y<a fl <27rJ m D l7 -<a D <y in D 2l (27) 
where D 1 , D 2 C <9f2 are disjoint sets with positive Hausdorff measure. 



5.3 The full problem in the elastic regime 

Now we investigate a case complementary to the one previously studied where 
the material is so hard that no plastic flow occurs. We assume again 70 = 
kq = 0. Computing the first variation of Ep given by f[T§j) w.r.t. 7, since 

f +1, if 7 >0, 
d sub | 7 | = < -1, if 7 < 0, 

{ [-1,+1], if 7 = 0, 

we obtain that no plastic flow occurs (i.e. 7 = 70) if 

— 2((/i c — A— n) cos a+A+/i) sina — /3(/i+(A+/i) sin 2 a+/j c cos 2 a) E [—ay, cry]- 

(28) 

A sufficient condition independent of a for (128]) is 

2 max{2(A + y) - /i c , ji c } + \/3(t) \ (A + 2/i + y c ) < a Y . (29) 
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For fixed A, fi and /i c , Condition (I2"9j) is satisfied if the yield stress cry of 
the deformed material is large enough and then, (T2^|) turns into a smallness 
condition on |/3(t)|. 

The rotations R c are determined as solutions of the minimisation problem 

Ep(a) := J 2/i 2 |Va| 2 + 2(A+/i)(l-cosa) 2 + 2/i c sin 2 a 



+2/3(t)((yU c — A— ji) cosa + A+/i) sina 

H (yU + (A + /i) sin 2 a+/i c cos 2 a) dx — > min 

2 

subject to (1201). 

When computing the Euler-Lagrange equation, we thus need to find a G 
L 2 (ft; [0,27r)) which solves 



\2. 

sin a 



— yU 2 Aa; + A+yU + yfi c — A — /it) (1 — Jcosa 

+ ^(A+yu) cosa + (yU c — A — yu)( cos 2 a — sin 2 a) j = in f2. 

(By standard elliptic regularity theory, then a G H 2 (fl; [0, 27r)) ). 

This time it remains to find stationary solutions to the Allen-Cahn equa- 
tion ([MD with 

Jp{a) := — I A+/i + p(t) sin cm cos cH (A+fijsma 

V 2/2 

H c -\-H( (3(t) 2 \ . 2 
H (1 — Ism a. (30) 



For = 0, Jp is identical to ( )25|) and the partitioning of Q into subsets 
with a ~ and a ~ 7r takes place if (126|) holds as discussed in the previous 
subsection. When (3^0, the discussion is more elaborate as Jp depends 
on 2 independent parameters and a. Since Jp(a) is particularly simple for 
(3 = ±2, let T > be a stop time and (3(t) be any continuous curve with 
(3(0) = -2, f3(T) = +2 that fulfils Q22) uniformly in < t < T. 

Fig. [2] displays a typical plot of the two distinct local minima mi, m2 of 
Jp for /9 G [—2, 2] computed with a one-dimensional Newton method, here 
for the special choice /i c = 6fi, A = fi which satisfies f f26|) . 
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Figure 2: The two local minima mi((3), m 2 (/3) of Jg as a function of /3 G 
[—2, 2] for /i c = A = /i. 

Like the condition f|2T|) . we can derive sufficient conditions on for the 
formation of boundary layers. A strong time-independent condition is 

(o<a D < 71^(2)) or (m 2 (-2) < ai -,<27r) in D h 

Mi(-2) < a D < M 2 (2) in D 2 . (31) 

where Di, D 2 C <9f2 are disjoint sets with positive Hausdorff measure; 
Mi(t) < M 2 (t) are the local maxima of Jp, and it was used that M\(t) 
is strictly increasing, M 2 (t) strictly decreasing in time. 



6 Discussion and concluding remarks 

In this article we computed analytically solutions of a 2D shear problem 
within the framework of finite-strain Cosserat plasticity. 

As main result of our investigations, we considered two complementary 
cases of shear in 2D, one of an ultra-soft material where plasticity occurs 
during the entire deformation process, and another case of a hard material 
where the applied loads are not large enough to initiate plastic flow. In both 
examples it was shown that the micro-rotations R e are stationary solutions of 
the Allen-Calm equation ([23]) which may result in a partitioning of the mate- 
rial into cells due to the occurrence of transition layers in R e . If y/JI^ is small 
(compared to |J1|), the transition layers are steep and R e is approximately 



13 



constant in each cell with a value determined by one of the local minimisers 
of a functional J that depends on the parameters A + /x, (fi c — A — /z)/2 and 
the applied deformations f3(t), cf. Eqn. fl30l) . 

The simple example in Subsection 15. II demonstrates that this partitioning 
is restricted to a certain parameter range of 1x2 and does not always occur. In 
addition, like the celebrated Cauchy-Born rule, the values of R e at dfl deter- 
mine the behaviour in the interior, especially whether and where transition 
layers occur. 

Yet, in many cases the analysis suggests a subdivision mechanism of Q 
caused by deforming the solid that may lead to an explanation why grains 
and subgrains form. Along this line, it is essential to work out the precise 
relation between local orientation R e and orientation of the subgrain. A first 
investigation in 3D using finite-element computations is done in [8]. In this 
work, it was assumed that both coincide. 

Since the analysed model relies on simplifying assumptions, it is natu- 
ral to ask whether the proposed mechanism carries over to the full three- 
dimensional setup, to general deformations, and to more realistic dislocation 
models. Another very important question left to future work is the study of 
e-minimisers, characterised by 

E(ip £ , R eE , F Pe , K e ) <m + e, 

where m is the absolute minimal value of E and e > a small number. 
This notion goes along with the insight that most of the time, the solutions 
observed in experiments do not reach a global minimum. If an almost optimal 
micro-rotation R 0£ solves the time-dependent Allen-Cahn equation (1241) for 
large but finite time t, this may pave the way for quantitative studies since 
the coarsening laws of the Allen-Cahn equation are well known. 
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